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The previous theory is extended to the case of an asymmetrical potential of the main body especially in 
the case of small eccentricity by assuming that J n , the coefficient of the »th-order zonal harmonics of the 
potential, is of «th order of magnitude and the eccentricity is of the first order of magnitude. Here there is 
a peculiar kind of libration which never occurs in the case of moderate eccentricity. This peculiar kind of 
libration splits into two kinds of libration, depending either on the antis ymmetrical term prevailing case 
or on the symmetrical term prevailing case , which is a continuous transformation of the type described in the 
previous theory. Numerical test discloses that for the earth the former peculiar kind of libration occurs. Also 
it is shown that the fifth coefficient is comparatively large and plays an important role in the asymmetrical 
theory for the case of small eccentricity. The present paper shows that the antisymmetrical terms cannot 
be neglected for the earth in the case of small eccentricity. ^ . /* 
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I. INTRODUCTION 

TN the previous paper by the same author (1962) the 
motion of a dose satellite in the vicinity of the 
so-called critical inclination was studied. However, in 
that paper the author omitted the effect of the anti- 
symmetrical terms in the potential of the main bodies. 
The theory including these terms was given, for 
example, by Kozai (1961). However, in spite of his 
comment on the case of small eccentricity, it does not 
seem very extensive. 

Therefore, the author has decided to develop a theory 
to avoid the difficulty connected with small eccentricity. 
A preliminary consideration shows that this difficulty 
arises only in the vicinity of e=0, and that there is no 
trajectory in general extending from the vicinity of 

* This work was performed under a National Academy of 
Sciences postdoctoral resident research associateship program con- 
nected with The National Aeronautics and Space Administration. 

f Present address: Tokyo Astronomical Observatory, Mitaka, 
Tokyo, Japan. 


e=0 to a value of the order of unity. Keeping this 
in mind, the author has expanded the Hamiltonian 
into a power series in e, which is assumed to be of the 
first order of magnitude, assuming that 7 n , the co- 
efficient of zonal harmonics of nth order, is of the nth 
order of magnitude, and that the deviation given by 
the following formula is assumed to be of the second 
order : 

a=l-SH 2 (fiao)-\ (1.1) 

where H, a constant, is the projected angular momentum 
to the equatorial plane, and a Q the mean semimajor axis. 

The necessary terms up to the sixth order of magni- 
tude in this respect are picked up from the Hamiltonian 
F*, which is given by removing the so-called periodic 
perturbation terms. The terms lower than the sixth 
order are all constant, therefore the sixth-order terms 
are considered as the leading terms in our theory. 

The author found that the equation of motion can 
be reduced to the same type as the ones given by 
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Andoyer (1903) : 

dh dF 

dr dh 

where 

F=-0k+(y+l3')h* 

with a different restriction on the sign of ft than that 
given by Andoyer. 

He assumed a restriction such as ft>0, which is 
always confirmed in his theory connected with the 
librational problem of asteroids. In our theory, how- 
ever, this restriction should be removed; namely, in 
some cases we have ft<0 as well as ft>0. Especially 
in the case of the earth, ft is negative. Therefore, some 
change from his theory must take place. This slight 
change might easily be overlooked, but as seen later in 
this paper, some alterations are required afterwards, 
if we want to have real expressions for the solutions. 

The main purpose of this paper is to develop the 
theory with this difference in mind, and exclusively to 
give the case for ft <t) . Section II presents a prelimin- 
ary discussion on obtaining the equations of motion in 
Andover’s form. Section III is the same as that of his 
theory; however, in order to avoid some confusion the 
present author rewrites the results in a very compact 
form. In Sec. IV some changes from his own appear in 
separating the several cases connected with the relation 
between the quantities y* = y/ft and ft' = ft/ ft. 
Section V gives a classification under which the solu- 
tions connected with the quantity n^ — n/ft, where u is 
the energy constant of the system, should be written 
down separately. Section VI gives the whole expressions 
of the solutions in the real number representation. 
Section VII is an Appendix which gives the char- 
acteristics, which means the trajectory within the plane 
of e cos g and e sing without any attention to the relation 
with the time, the independent variable. Several per- 
tinent numerical results are given. Section VIII gives 
some discussions related with the convergency, the 
relation to the previous theory and so forth. Finally 
Sec. IX gives the conclusion. 

In order to test the assumption imposed on J n , the 
quantities for the earth are listed below — not only in 
actual values but also in units of proper powers of 
aojd»— 0.0528, which is assumed to be of the first order 
of magnitude : 

J 2 = 1082.36X 10- 6 = 0.3882 0 2/ W, 

/ 3 = - 2.566X 10“ 6 = —0.0174 £ 3 % 0 3 , 
y 4 =_2.14 x 10- 6 = —0.275 ft*a 0 \ 


dk dF 

dr dk 


( 1 . 2 ) 


f {y-ftW+{h 2 +ft)\ (1.3) 


J 5 = -0.063X 10~ 6 = -0.154 

Roughly speaking, these numerical values show that 
/ 2 , J 4, and 7 5 play approximately the same order of 


role in the vicinity of 6=0, say 6=0.10 or 0.05; on the 
other hand, / 3 does not take any important role here. 

II. EQUATIONS OF MOTION 

The potential, under which influence a negligible 
small mass particle moves, is assumed in the following 
form : 


Mr 1 , v 

r =- 1 -T, Jny—J Pnisind) , ( 2 . 1 ) 


where v=Mk 2 , k being the Gaussian constant and M 
the mass of the main body, r is the distance of the 
particle from the center of the main body, 5 is the 
declination of the particle, /„ represent numerical 
constants which characterize the spheroidal potential; 
in this paper the summation £ extends from 2 to 5, as 
seen later. 

Using a result of the so-called secular and long 
periodic parts in the original Hamiltonian such as 
given by Kozai (1962a), we may pick up only the 
following necessary parts provided that /„ is assumed 
to be of the nth order of magnitude as well as the 
eccentricity e to be of the first order: 

3 ft J 2 21 ftj a 3 ftJ 2 2 

F* = (e 2 — a) 2 ~\ e 2 e 2 cos2 g 

16 TJ 40 D° 40 L 10 


3 ft J 4 


40 V 


-e 2 cos2 g- 


3 ftJz 


4(51) l 8 


-e(e 2 — a) sing 


3 ftJ 2 Jz 
4(5*) L 12 


-e smg- 


9 ftJ, 
10(5*) L 12 


-e sing, (2.2) 


where a= 1— 5H 2 (nao)~ 1 = a. constant which is assumed 
to be of the second order, and the other notations 
correspond to Delaunay’s. In this expression, the terms 
of beyond the sixth order of magnitude are neglected. 
As the Hamiltonian has only sixth order of magnitude 
and nothing else, neglecting the higher orders, we may 
take £ and rj as the canonical variables, 

Z*6 sing = £, Z,*6 cosg = 77, 

in the canonical equations of motion : 

d( dF * d v dF * 

-=— • (2.3) 

dt dr) dt 


For mathematical simplicity, if the canonical 
variables and the independent variable are changed : 


k — e cosg, hi=e sing, 
3 ftJ 2 3 ftJ 2 
16 V 16 a 0 712 ’ 


(2.4) 
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then the equations of motion will be 

dh\ dP \ dk dF \ 

dr dk dr dh\ 


(2.5) 


where by neglecting an unnecessary constant in F *, 


_ 16 L 6 

Ft — (F* — const) 

3 fjdJ 2 

= — pih+ (7i+jSiOAi 2 + (yi—pi)k 2 

+/3 x 'W+£ 2 )Ai+(Ai 2 +£ 2 ) 2 , (2.6) 


with 


4/ 3 4 Jz 24 / 5 


j3i=+cr 


5*/ 2 a 0 5*a 0 3 5(5*) / 2 a 0 3 
14 J A 


7i = - 


-2a, 


/3i" = 


5 / 2 a 0 2 

2/ /A 1 
1 '=- (/*+-} — 
5\ /*/ a 0 2 

4 J z 


5* J 2&0 


(2.7) 


In order to remove the third-degree terms in hi and k , 
Ax is changed to h given by 

A-Ax+W'. (2.8) 


Further, if the equations are transformed in terms of 
p and q , 


p=h+ik, q—h—ik , 

(2.9) 

the equations can be written as follows : 


dp dF dq dP 


— =-27— , -=+27—, 

(2.10) 


dr dq dr dp 


and 





The form of the Hamiltonian (2.11) shows that it is 
the quadratic in p or q, respectively; this was the 
technique which Andoyer (1903) used for the problem 
associated with the libration near the commensurability 
between the mean motion of asteroids and that of the 
disturbing body, Jupiter. Therefore, we can follow 
after his development in order to solve the equation of 
motion. Nevertheless, there is a slight difference be- 
tween his case and ours concerning the sign of some 
coefficient. He assumed that the coefficient /3' is always 
positive; but in our case can take a negative value 
as well as a positive one depending on the interrelation 
of the magnitude of J n * s. Especially for the earth it is 
negative, as seen later. Therefore, we are restricted to 
take the negative /3' case here unless otherwise men- 
tioned, because for the positive case Andoyer’s results 
are available. 


III. GENERAL EXPRESSION OF SOLUTION 

The general form of the solutions does not change 
from Andoyer’s results. But, in order to avoid some 
confusion about the real/complex criterion, we shall 
make a slight change in the notation. The calculations 
are omitted and only the results obtained are given here. 

The solutions of Eqs. (2.10) are given in the following 
form: 


1 

p= [f (r+X+|p— ( r +^ — hp~hx) 

7(a 0 ')* 

— if(p+x) — U(p~ x)], 

1 

?= [f(r+X+£p+£x)— f(r+*“ip+$x) 

7(a 0 ')* 

-it(p+x)— 2.f(p— x)]> 


with 

F=—?&(P+q)+ypq+hP , (iP—q 2 )+iP<f> ( 2 - 11 ) 

where 


„ _ 2 7 3 2 18 7 3 2 16 7 3 2 7 4 

F-Fi= a h 

5 7 2 W 25 7 2 U 0' 1 25 7 2 3 ao 4 


+ 


24 J 3/6 3 7s 4 


+ 


25 J 2 W 25 J 2 W 


-= const, 


/ 4 8 7 4 2 7 3 2 \ 4 7 3 24 Jf, 

0 =(— + ) , 

\ 5 5 7 2 2 5 /•// 5» a<? 5(5»)7*«o* 


( 2 - 12 ) 


with 


p(p) = ~ (ai 2 +ao'a 2 )(a 0 ')- 1 (real), 

fp' (p)—i(ao' 2 a3+3a 3 'aiai+ 2a v 3 ) (a 0 ') _l (imaginary), 

^(x)=8?(p)+ii(ao') J |3-V'(p) ( real ), (3.2) 

p'(x) = -iS'jS-^floOV'W (imaginary), 

where 


= — 8/3', d\— 2/3, 

a 2 =K 8 «+ 2 ^' 2 - 27 2 ), a z =p(y-(3'), (3.3) 

a 4 =8wj3'— /S 2 ; 


14 7 4 4 7/ 

7 = 2aH - , 

^ J 2 a 0 - 5 7 2 W 


u is the energy constant such that 

F+u= 0 . 


(3.4) 
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#>(z) is the Weierstrassian elliptic function with the 
parameters 

#2= — at/04—4aia3+3a2 2 , (3.5) 

ai 2 a4+uoW— ai, 


and so forth, then 


A* = 1024<2* 2 y.+p^u* 1 


such that 


Lfp'(z)y=-ip s (z)-g-z&(z)-g. h 


where p'(z) denotes the derivative respective to z, 
f(z) is the associated zeta function defined by 


fW 'H 


In the expression of solution (3.1) X is an integration 
constant, which with the energy constant u forms a 
system of the arbitrary constants in the solution of the 
equations of motion, u is of course a real value; X is, 
however, not necessarily real. X should be taken such 
that the expression (3.1) could give real values of h 
and k , ~ 


h = i(p+q), k = — (p~q). 

2 i 


IV. DISCRIMINATION AMONG SEVERAL CASES 

At the first step it is necessary to know the sign of 
the discriminant A, 




For brevity, let 


+ it(t*+^*0L(7*+^+0 3 +9]w+ 

- ^[ (7 * + /5* ,)3+ t]1' (4.1") 

If the subscript * is omitted in this expression, then 
every term remains the same as in (4.1') except for p 
which is replaced by unity ; besides, as is easily seen, 
P is of the third order of magnitude, and we may choose 

h*=h/fa k*=k/pK (4.4') 

As a result, all quantities are measured in units of 
proper powers of £*, which is assumed to be of the first 
order of magnitude, provided that Therefore 

without any restriction we may take p=l hereafter 
unless otherwise stated. This is an application of the 
nondimensional analysis. 

In any way, there are two cases for A=0: 


If we consider (4.5) as an equation for w, then we 
know the value w 0 which satisfies this equation: 

«o=i (P'-yT+W- 1 . (4.6) 


(ii) /(«) = u 3 2 (y+P ')«H [(t+i8') 3 +9]m 

16 


P= P(x)=i( 4 «— yt+bp’y— 5 / 9 ' 2 ), 

and (4.2) 

Q= 16w/3'-4 P'iP'-y)*- 2p, 

which are also real values even for our case, then 
fl -12P*-4( 7 -2 (3')Q, 

and (4.3) 

gi =-&P s +Q*+My~2p')PQ. 

It follows accordingly that 
A = 1024@ 2 j m 3 — § (y+i8') 2 w 2 


ir 27-] 

(7+i3') 3 =0. (4.7 

64L 4 J 


Here two cases again are divided according to the sign 
of the discriminant A' of f(u) itself: 


1 f 27 

A'= (7 +PY+~ 

1281 8 


namelv, if 


7 +£'>“§, 


+-&{y+P')\_{y+P') i +W]u 


—hP (y+P7+~P^. (4.1') 


For brevitv, let 


«*=«/ P* 9 8*=P'/ p* 9 

A*=A/0 8 , 

Q*=Q/P 2 , p+=p/P\ 


then among the three roots u h w 2 , w 3 of Eq. (4.7) in 
terms of «, two of them have complex values, while 
the other, say u h remains real. On the other hand, if 

7 +P'<~1 (4.10) 

then the three roots are all real, 

U-Z<U2<Ui. 

Now, it is necessary to know in what sequence the 
roots «t(i=0,l,2,3) appear; in other words, to know 
the relation between and » t *(i= 1,2,3)(A'>0) or 
Mi(A'<0). In the case when (4.9) occurs we can easily 
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determine it in the following way: first of all we have 
/ («)^0, if and only if wig u \ ; (A) 

however, since 

9 

/(« 0 ) = -(7 2 +^'- 1 )Po 2 , (4.11) 

64 

where Po denotes the value of P when m — u 0) we have, 
if /3'>0, the right-hand side of (4.11) is always positive, 
but if 0'<O, this is not always true. Consequently, we 
have a new discrimination, which never occurs in the 
case when p> 0: 

if y 2 +%&'~ l >0 then u^ui. (B) 

On the other hand, in the case when (4.10) is satisfied, 
we may only say that 

if 7 2 +^ ,_1 <0, then « 2 <w 0 <«i, 

or u 0 <uz; (C) 

and if 7 2 +i 0 / ” 1 >O, then Mi<«o, 

or i<3<Wo<w 2 . 

The detailed criterion will be given in a later part. 

At the second step it is necessary to divide the cases 
according to the sign of (g 3 ),- (*= 0,1, 2, 3), where (g 3 )» 
denotes the value of g 3 when u—Ui. This is important 
because, for example, if (g3)t<0 then the energy con- 
stant w; corresponds to the unstable equilibrium 
point (s). It is easily confirmed for (g 3 )o, since 

512 

(g 3 )o= -8Po 3 = — ^ 3 (7-^+i^- 2 ) 3 ; (4.12) 
27 


when /3'<0, 

if 7-0'+i/3'- 2 <O, then (g,) 0 < 0. (D) 

In order to know the relation between u 0 and 

(i = 1,2,3) in more detail as well as the sign of (g 3 ) t - it is 
convenient to divide the 7 vs fi f plane into several 
parts so that in each of them the situation never 
changes. For this purpose, it is sufficient to draw lines 
where the situation would change. The following four 
lines serve this purpose : 

7 +£'=-§, (4.13) 

7+(-2/?')-* = 0, (4.14) 

7-(-2/3')-*=0, (4.14') 

and 

7-/3 , +i(/3 / )“ 2 =0. (4.15) 

In effect, Uo—Ui occurs only on the lines (4.14) or 
(4.14'); U\ = Ui or w 2 = w 3 only on the line (4.13); and 
(gs) o=0 only on (4.15), (g 3 ), = 0 (i^0) occurs on 
either (4.14), (4.14') or (4.13), as is easily confirmed. 
Before dividing the plane by the four lines given 


above, for the sake of simplicity let us consider the 
division by separating the negative ft' half-plane into 
three parts : 

Case I. 0 '<~ 2 ; 

in this case we have 

-^-f> + (-2d')-i>-(-2/?')^>/3' 

-l(m (4.16) 

Case II. -2 <d'<-f, 

(-2 £')-*> — /3'— f > - (-20')“*>0' 

-urn*- (4.i7) 

Case III. — \ <^,^'<0, 

(— 2/3')”i> — — (— 2/3')~ i >/?' 

-WY 1 - (4.18) 

The reason why we make a discrimination between 
case II and case III is that 

if 

then 

- (-2/?')-*=tf'-itf'‘ 2 = - 1, 

so that a domain, — — § > 7> — (— 2#') - * for example, 
is cut out at £?'= — \ and there is no continuous route 
connecting the regions corresponding to case II and 
case III without meeting any one of the above four 
lines. 

Now, we subdivide the regions with respect to the 
four lines stated above; for example, case Ii means 
0'< — 2 and 7 >— £'—f, case I 2 ft'< — 2 and — 

—%>y(—2p'Y and so forth. This subdivision is not 
essential (there is, nevertheless, a practical advantage, 
because 0', or speaking more precisely = de- 
pends only on the coefficients of the zonal harmonics 
for the potential and is totally independent of the 
initial condition) as a whole because cases Ii, Hi, and 
IIIi, for example, have no differences in the sense of 
continuous deformation. Therefore, it is useful to 
assemble some of the too subdivided cases into one 
case, then finally we have the foUlowing division : 

Case 1 which involves cases Ii, III, IIIi ; 

Case 2' which involves cases I 2 ; 

Case 2" which involves cases II 2 , III 2 ; 

Case 3' which involves cases I 3 , II3 ; 

Case 3" which involves cases III 3 ; 

Case 4' which involves cases I 4 , II4 ; 

Case 4" which involves cases III4; 

Case 5 which involves cases I 5 , II 5 , III5. 

Figure 1 shows the division in the fi' vs y plane. 

Here for the sake of comparison, the case where jd'>0 
(Andover’s subcases a', a ", b\ and 6") is added. 

We evaluate tu and (g 3 )» in each case either by 
giving extreme values for ft' and 7 or by giving special 

i 
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y * 



Fig. 1. Division in the /3'+ vs 7 * plane. and 7 * = 7 //3 l . 


numerical values as given in Sec. VII. In any case one 
set of values for and 7 is enough to determine the 
sign of (g 3 )t and the relation of 1/0 and Vi ( i= 1,2,3). 

As an example, case 3' is shown where extreme values 
are given. 

Case 3' 0'=-e-*, *->+0, 

512 

($3)0= H 

27 

2 


and 


512 

(£ 3 ) 1 , 2 = H e” 2 * 

27 



where unnecessary higher-order terms are omitted. 
Therefore 

«3<tt2<S/0<t/l, 

and 


(gs)o>0, (g 3 )i>0, (g 3 ) 2 >0, and (g 3 ) 3 <0. 

Table I shows the results. From this table, combined 
with (4.7), it is easy to find the sign of A(w) [e.g., in 
the case 3', if wi<w then A(«)>0, if w 0 <w<t/i or 
w 2 <w<w 0 then A (ft) <0, and so forth]. 


V. VALUES OF p AND % 

For the purpose of giving a real expression of the 
solutions, it is also necessary to have the range in 
which p or x falls. As a first step, a comparison is made 
here between p(p) or p(x) and the parameters of p 
function, namely e h e 2, and e 3 . Let 

4^(z)-g2^(3)-g3 

= 4[ j?(z) - ei][> (z) - e 2 ][ j? (z) - e 3 ], (5.1) 

then it is well known that if the discriminant A>0, 
the three parameters e 1} e 2 , and e 3 are all real; contrary 
to this, if A<0, then only one of them, say e h is real. 
Since 

&>'(p) 2 = 4[>(p) ~ «i][>(p) - ^DEXp) - es], (5.2) 
%>{p)= — (ai 2 +ai 2 a 2 )flo' is real, (3.2) 

and 

g?'(p)“f(ao 2 a 3 +3^o / tii«2+2ai 3 )(^o / )' _! is imaginary, 


Table I. Roots of A(m) =0 and sign of g* at respective roots. 


Case 

Ui 

(gs)o 

(ft) 1 

(ft). 

(ft). 

Case 1 

Uy<Uq 

+ 

+ 



Case T 

Ut<U 2 <Ui<U 0 

+ 

+ 

+ 

— 

Case 2" 

Uq<Ui 

+ 

+ 



Case 3' 

U1<U2<U0<U1 

+ 

+ 

+ 

— 

Case 3" 

U0<U1<U2<U1 

+ 

+ 

— 

+ 

Case 4' 

U1<U0<U2<UY 

+ 

+ 

+ 

— 

Case 4" 

U1<U0<U2<UY 

+ 

+ 

— 

+ 

Case 5 

UZ<U0<U2<U1 


+ 

+ 

+ 


it follows that, when A>0: 

p(p)<e z or e 2 <p{p)<e h (5.3) 

where the three arguments are denoted such that 

e 3 <e 2 <£i. 

Similarly, with 

& f (x)=P(p)+bM^ 1 p / (p) = P is real 
p'(x)=-P'0~ 1 M*p'(p) = iQ is imaginary, (4.2) 
when A>0: 

p(x)<es or e 2 <p(x)<ei . (5.4) 

In the case when A<0, it follows necessarily that 

and (5.5) 

p(x)<ei- 

Besides, from (4.2) it follows that 

&(fi)>&(x), if Q>°, 

that is to say (/3'<0), (E) 

W:<Wo. 
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Table II. The values of P= g?(x) and p(p) (p'<0). 


Case 

Range of u 

Values of P and p 

Sign of 

and 
*p'(p) /* 

Class* 

Case 1 

u Q <u 

P<e%<e 2 <P<e x 

— 

01 


U = Uq 

P = et=e2= p 

0 



Ui<U<Uo 

P<e i<e 2 < p<e x 

+ 

02 


U<Ui 

P<»<e 1 

+ 

12 

Case 2' 

Uq<U 

P<ez<e 2 <P<e x 

— 

01 


U = Uq 

p = e z = e2 = P 

0 



Ui<U<Uo 

P<ei<e 2 < p<e x 

+ 

02 


Uz<U<Ui 

P<V<e , 


12 


Ui<U<U2 

P< p<et<e2<e x h 

+ 

22 


«<W| 

P< v<ei b 

+ 

12 

Case 2" 

Ui<U 

$><ei<e2<P<ei 

— 

01 


Uo<U<Ui 

V<P<e 1 

— 

11 


U — Uq 

p=et=e 2 =P<ei 

0 



U<Uq 

P<9<e ! 

+ 

12 

Case 3' 

Ui<U 

$<ei<e 2 <P<e 1 

— 

01 


Uq<CuKui 

9<P<e 1 

— 

11 


U = Uq 

p=e%=e 2 = P 

0 



U2<U<Uq 

P<»<e 1 

+ 

12 


Uz<U<U 2 

P<»<ei<e 2 <e l h 

+ 

22 


U<Ui 

P<V<e i b 

+ 

12 

Case 3" 

Ui<U 

p<e z <e 2 <P<ei 

— 

01 


U2<U<Ui 

&<P<e i c 

— 

11 


U%<U<U2 

&<P<e z<e 2 <ei c 

— 

21 


Uq<U<Uz 

$><P<e 1 

— 

11 


U — Uq 

p = e 2 =et=P<ei 

0 



U<Uq 

P<$<e 1 

+ 

12 

Case 4' 

Ui<u 

%><ei<e 2 <P<e x 

— 

01 


U 2 <U<U X 

p<P<e 1 

— 

11 


Uq<U<Uz 

ei<e 2 < %><P<e 1 

— 

23 


U — Uq 

P = ei=e 2 =P 

0 



Uz<U<Uo 

P< p<ei<e 2 <ei h 

+ 

22 


u<u x 

P<V<e i b 

+ 

12 

Case 4" 

Ui<U 

p<ei<e 2 <P<e x 

— 

01 


U 2 <U<Ui 

p<P<e i° 

— 

11 


Uq<U<U 2 

V<P<e%<e 2<eP 

— 

21 


U~Uq 

p=P=e t =e 2 

0 



Uz<U<Uq 

e%<et <P< g><«i 

+ 

24 


U<Ui 

P< p<e x 

+ 

12 

Case 5 

Ui<U 

p<et<e2<P<e x 

— 

01 


U 2 <U<Ui 

V<P<e 1 

— 

11 


Uq<U<Uz 

e%<ei< $><P<e 1 

— 

23 


U=Uq 

e 2 =P= p = ei 

0 



Ui<U<Uq 

ei<et<P< p<e x 

+ 

24 


U<U Z 

P<P<e 1 

+ 

12 


• See Table III. 

b ei in the lower line is the continuation of es in the upper line as a result 
of the dropping of ei and e 2 in the upper line into imaginary values. 

0 ei in the upper line is the continuation of e\ in the lower line as a result 
of the dropping of ei and et in the lower line into imaginary values. 


The above is simply obtained. The next is to decide 
whether p(x) and p(p) drop either in the region 
between e\ and e 2 or the region smaller than e 3 . This 
decision can be given by knowing the behavior near 
where generally the situation is changed, and 
by knowing the behavior in the extreme case such as 
u — > + 00 . The detailed calculations are omitted here, 
and only the results thus obtained are listed in Table II. 
In this table, for brevity, we let : 

p=p(x) and p=p(fi). (5.6) 

Besides, P'(x)/i and p’(p)/i are determined by 

&’(x)/i=Q= 16/3' (w— Mo), 


which give us the sign of both values, provided that 
(go')*, or more rigorously (ao')**= (floOV0*> is non- 
negative without any restriction. In effect, since in the 
form of the solutions h and k, or more rigorously 
h*=h/0^ and k*=k/0*, have (ao') -i , or more rigorously 
(goO *“* = /?*/ 2 ( — 20') as a factor, the change of sign 
of (a o 0* as a whole does not produce any real difference 
in the solutions at all. It is to be noted that when 0<O 
this provision means (ao') _i itself, which does not 
mean a nondimensional quantity should be 

taken as negative. By this method we may ignore the 
difference of sign of 0 but may unify the procedure as 
a whole. This is justified directly from the fact that in 
the equations of motion the sign of 0 can be changed 
without any significant alternation except for the change 
of signs of h and k for the same time. 

The next step is to determine the range in which p 
or X will fall. For this purpose it is convenient to 
classify the various cases which are associated also 
with the value of u. This was done in the last column 
of Table II. The specification of each class is given in 
Table III. Thus we can determine the range of p and x, 

if we remember that : 

when A <0, if 0O<a)i", then 

“°° <$>(iv)<ei, and g>'(iv)/i< 0, 
but if coi" < v < 2a>i", then ^ 
ei> — °o, and fp'(iv)/i> 0, 

where c oi'i is the purely imaginary semiperiod given by 


COi 


tt _ = 



dp 

[4 {ei-p) (e 2 - p) (ez-pytf 


(5.9) 


when e x is real. 

On the other hand, when A>0, 

if 0<p<o/', then 

oi Jp (uJ “k M*) ^ £2 J p (iv f ~J~ iv)/ ^ 0, 

— <X> <fp(iv)<e 3 , and p'(iv)/i<0, 
but if «" < v < 2a)", then 
e 2 <p(o> f -\-iv)<e 1 , gp , (o/+ii')/i<0, 

e z > p(iv)> — *>, and p r (iv)/i> 0, 


where a/ is the real semiperiod and a)i"i is the purely 
imaginary semiperiod given by {ez<e 2 <e^) 


j 


dp 

[4(/>— ei) (p—e2)(p—e 3 )2 i 


(5.10) 


/ tz 

-00 


dp 


[4(«i -p)(ei-p)(e 3 -p)y 


and 


p'(p)/i= - 16 (ao')~*(«-wo), 


(5.7) 

(5.8) 


and 


(5.11) 
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Table 111. Specification of class and range of p and x- 


Class 

Specification 

Sign of v'{p)/i 
and of *'&)/* 


Range of p and x 

01 

p<e3<ei<P <ci 

_ 

P = 
ip, 

X ~ 

w'-Hx' 

ix 

0 < p '< m "< x '< 2 u " 

0<p'<w''<x'<2&>" 

02 

P<e 3 <e 2 < &<ei 

+ 

0)' + lp', 

11 

&<P<e \ 

— 

ip', 

ix 

0<p'<x'<W' 

wi"<p / <X / <2coi' / 

12 

P< $<e x 

4- 

ip', 

ix 

21 

y><P<e 3 <€ 2 <ei 

— 

ip', 

ix! 

0<p'<x'<<*>" 
(*>"<p' <x 

22 

P< p<ez<e 2 <ei 

+ 

i P , 

ix' 

23 

e 3 <e 2 < %><P<e i 

— 

w'4 -ip , 

«'+*x' 

w'-Hx' 

(*>" <p' <X* <2o>" 
0<p' <x'<o>" 

24 

e 3 <e 2 <P< £><ei 

+ 

Note: Class 0 and Class 2 are corresponding to A >0. On the other hand. Class 1 is corresponding to A <0. 


VI. 

REAL EXPRESSION OF SOLUTIONS 

canceled out 

by choosing 

the proper value for the 


As seen in the preceding section, the values p and x 
are uniquely determined within the parallelogram, 
which is (2a/,2w "i) for A>0 or (2co/,o;/+a j/'i) for 
A<0; however, the parameters which are included in 
the expression of solutions are not p and x, respectively, 
but fp and ix* These values can not be uniquely 
determined in the parallelogram but correspond to the 
four values, respectively. In effect, v is uniquely 
determined by: p(v)= given constant, has given 

sign, then 

ip+to', 4v+w"?‘, and i^+a /+«"*, for A>0, 


or 


and 


ij'+te/, ip+i^Z+co/'i), 


ff i f for A <0, 

are all required values as well as \v. 

This phenomenon, however, does not provide any 
essential difficulty. In effect, the change from \v to 
i^+a / produces only a half-period advance, and this is 


additive constant X. On the other hand, the change 
from \v to or to produces 

another series of solutions, which satisfy the equations 
of motion as well as the given energy constant v. 
Reviewing the results in advance, class 2 (which in- 
cludes classes 21-24) has two series of real solutions; 
contrary to this, class 1 has only one series of real 
solutions while the other corresponds to imaginary 
solutions, and class 0 has no series of real solutions. 
This change can be also carried out by changing X to 
X'+a/'i or X'+i(wZ+<o"Z) without any alteration of 
^p and ix, respectively, as was done by Andoyer, since 
the solutions are periodic qua function of X, |p or ix, 
respectively, and the argument is in the form of Xdt|p 
zb^X. Accordingly, we may fix the values for ip and ix. 
And from 

§(p+x) = Xi, 

and (6.1) 

Up— x) = X2, 

we may construct Table IV. 

Thus we have obtained the real expressions of the 
solutions in each class. From (3.1) the expressions of 


Class 

01 

02 

11 

12 

21 

22 

23 


Table IV. Constants Xi and x>. 


Xl — 

i(p+x) 

X2 = 

i(p-x) 

— Jw'-h Xi% 

W' <x! <W" 

= ~w~\-x a 

1 

A 

x 

A 

o 

sw'+x/f, 

W'<x!<W' 

\i*>'-\-Xi'i 

, <X»'<0 

Xi'i, 

0<Xl'<«l" 

x!i, 

o 

V 

X 

V 
3 
1 

Xi'i, 

o 

V 

* 

V 

3 

1 

x/i, 

— W'<x 2 '<0 

Xi'i, 

0<Xi'<w" 

x!i, 

-iw"<X2'<0 

Xi'i, 

w"<Xi'<2a/' 

Xi'i, 

— W'<x!<0 

w/+x/*. 

w ;/ <x/< 2to" 

x!i, 

— W'<x!<0 

+ Xi'i, 

0<Xi'<p>" 

Xi'i, 

l 

a' 

X 

A 

o 


24 
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Table V. Constants included in solutions (6.8,9). Xi=Xi'*\ X 2 =X 2 '^ 


Class 

Specification 


Sl 

CQ 

A 

X," 

x*" 

xr 

/// 

x 2 

Hi 

co" 

0<x/< — , 
2 

OJj" 

<x/<0 

2 

-1 

-1 

+1 

0 

x/ 

2 

o>i" 

x*'+— 

2 

2xi' 

2 

0>l" 

2x 2 'H 

2 

11* 

0>l" 

<Xi'<o>i", 

2 

0>l" 

<X 2 '<0 

2 

-1 

-1 

-1 

+i 

0>l" 

Xi' 

2 

»i" 

Xs'H 

2 

2xi'— §«/' 

Wl" 

2x«'-f— — 
2 

12 x 

1 

A 

x 

A 

© 

0>l" 

<X 2 '<0 

2 

+ 1 

+ 1 

-1 

+\ 

«i" 

x/4 

2 

«l" 

X 2 '-h— 

2 

2x/H 

2 

0>l" 

2x 2 '4 

2 

12* 

o>i" 

— Wi' ' <Xi' < > 

2 

Wl" 

<x 2 '<0 

2 

+ 1 

+ 1 

+1 

0 

on" 

Xi'H 

2 

COi" 

X 2 '-i 

2 

2xi'4-fo>i" 

an" 

2x/H 

2 


Note: The subdivision here is caused by the fact that the restriction |Im(»')| <an" would not be satisfied especially for f(2xi) in the second subclass; 
accordingly before expanding into series we have changed the argument keeping in mind the periodicity of h and k qua function of xi- 


h and k are easily obtained : 

h=Up+q) 

1 

_ Cf ('J — hXH-Xi) — f (t+X— X x ) 

4^-20')* 

+r(r+X+X 2 )-f(r+X-X 2 )-f( 2X 1 )-J’(2X 1 )] (6.2) 




1 

= * [f (r+X+Xi)+f (r+X — Xi) 

4(-20')‘ 

— f (r+X+X 2 ) — f (r+X— X 2 )]. 

When A>0, the expression of f (v) is given by 

K)v 7 ri 2w oo Cj n me 

f W = — =F — +— L sin — (vTw"»), (6.3) 

C o' 2co' a/ n-1 1 — Q 2n Co' 


where the double signs should be taken according to 
the sign of the imaginary part of v : 

upper sign for 0 < Im(v) < 2co" 

or 

lower sign for — 2co"<Im(?') <0, (6.4) 

q=exp[— W'/w'], »)=?(&/), 

co' is the real semiperiod, and o o"i is the purely imaginary 
semiperiod; when A<0, 

t)xv x 

fW- =F i 

co/ 2c C\ 

2t Qi n mr 

H 51 sin — -F§coi"7) , (6.5) 

coi' l-cp 2n cox' 

where the double sign should be taken according to the 


Table VI. Constants included in solutions (6.10, 11). 


Class 


Specification 


Si 

A 

Xi" 

X 2 " 

Xi' 

// 

X 2 '" 

21 1 

Xi = 

Xl i, 

0<xi'<«i"1 

, ! 

+ 1 

+ 1 

0 

X,'-a>" 

X 2 '+«" 

2xi'- 

ai" 

2x 2 '+co" 

21 2 

X 2 ~ 

X«'i, 

02 

- — KX^O | 

+ 1 

+ 1 

0 

Xi' 

X*' 

2x/- 

a>" 

2x 2 '+a 5 " 

22 1 

[*- 

Xi'*, 

o>"<Xi'<2«") 

+ 1 

+ 1 


Xi'-oj" 

X 2 '+co" 

2x/- 

3o>", 

2x 2 '+oj 7 

22, 

[ X2= 

X 2 'a 

a; 

<x,'<0 

+1 

+i 

i 

4 

Xi'-2co" 

x 2 ' 

2xi'- 

3»", 

2x 2 '+cc" 

23i 

fx,= 

w'+Xi'i, 

d>"<X,'<2«"l 

n 

i (-1)" 

(-0” 

+i 

Xi'— w" 

x*'+»" 

2x/- 

■3w", 

2x 2 '+a>" 

23, 

l X,= 

x 2 'a 

<0 

<x 2 '<0 

2 J 

f (-i)» 

(-i) n 

1 

4 

Xi'-2co" 

X*' 

2xi'- 

■ 3oj", 

2x 2 '+ai" 

241 J 

rx-.= 

»'+Xi'«, 

0<Xi'<a>" ' 
// 

i (-i)" 

(-!)" 

0 

Xi'— oj" 

X 2 '+a>" 

2xi'— 

a>" 

2x 2 ' 4-05 " 

24, : 

l X!= 

x 2 'a 

05 

<X 2 '<0 

2 


(-1)” 

0 

Xi' 

x/ 

2xi'- 

-a?" 

2x 2 '+oi" 
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means the summation extending only over positive 
odd numbers, and only over positive even numbers. 
The coefficients are expressed by 

mrXi" mrX 2 ' 

GY = a' cosh h cosh , 

0)1 0)1 

mrXi" mrX 2 " 

C n ” = sinh bsinh , 

0)1 0 )\ 

mrXi" mrX 2 r 

S a — Si sinh sinh , 

0)1 0)1 

(6.9) 

mrXi" ' mrX 2 f 

Sn' = sinh sinh , 

0)1 0)1 

htXi" mrX 2 " 

C n0 ' = w cosh i cosh , 

COi / 0)1 

mrXi" < mrX 2 ff 

Cno"= — i sinh J sinh 

COi' COi' 


I 3.0 




The coefficients and the arguments involved in the 
above expressions are given for each class in Table V, 
respectively. 

A>0 (Class 0 has no real solution; on the other hand 
class 2 has two series of real solutions ; the existence of 
two subclasses corresponds to this phenomenon. The 
second subclasses are constructed by making X=X / +wY 
in the first subclasses, respectively. After being ex- 
panded into series, the unnecessary prime on X is 


omitted in the following.) 


1 7 r 

h= 

(-2/50*a/ 


r « c\ n ( nr ] 

X X) 1 C n cos — (r+X)+C’„o f +^4 

_n=l 1— q 2n l 0 )' \ 


1 7r * q" f mr 

k= X 1 S n sin— (r+X) 

(-2/9')* w' «=i 1— q 2n l to' 


where 


“ 7TCo""| 

q=exp , T=ysnKT 2 a<r' II2 L 

L o) f — 



( 6 . 10 ) 


Fig. 7. p ; *= — 1.0, 7 * = +2.0. (Case II i, or case 1 .) 
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The coefficients are expressed by 

mrXi' mrX 2 ' 

C n —c\ sinh hsinh , 

O)’ 0)' 

nnX i" nwX 2 " 

S n —si cosh cosh , (6. 1 1 ) 

w' 0)i 

mrXy!” nirX 2 "' 

C«o= sinh ^ sinh— . 

u' 0)' 

The coefficients and arguments involved in the above 
expressions are given in Table VI. These constants are 



chosen such that all the expansions should be con- 
vergent, or in other words, the restriction imposed on 
(6.3) is always satisfied. If this would not be satisfied, 
before being expanded into series the argument has 
been changed keeping in mind the periodicity of h and k 
qua function of Xu 

It seems necessary to add the dimensionality of co, 
Xi, X2, etc. They are all expressed in the following way: 
co* = a$*, and so forth. 

VII. CHARACTERISTICS 

Some numerical calculations connected with the 
characteristics are given here. The characteristics means 
the orbit without any regard to the independent variable 
r but only the plotted line in h and k space. In this 
problem it is easy to obtain such a characteristic. In 
effect, 

F+u = 0 (7.1) 



gives one ; or, in more detail 

(&+&)*+ ( 7 +&')&+ ( 7 -0')* 2 -^+« = O. (7.2) 

The same procedure of the nondimensional analysis 
taken in Sec. V gives the following form (omitting the 
subscript *) : 

(**+**)*+ (y+8'W+ (7-fl')**-*+«=0, (7.20 

from which we have simply 

h 2 — {^ 2 +i ( 7 00} 

±l- 2 p'k*+h-u+i(y-m. (7.3) 

Accordingly, bv putting 

U i= — 2/3'/z 2 +//+J(7— /?0 2 , 

and (7.4) 
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the following criteria are given : 

(i) if u> U i, then 0 solution for k ; 

(ii) i if «<U U A*+i( y-ft)>0, 

(ii)n and i*> t/ 2 , then 0 solution, 

(ii)i 2 or u<U 2 , then 2 solutions ; (H) 

(ii) 2 if u<U h A*+J( y-8')<0, 

(ii) 2 i and i/> Z7 2 , then 4 solutions, 

(ii ) 2 2 or u<U 2 , then 2 solutions. 

The critical values for u also come from these 
criteria as follows : 

/dUi\ 



Fig. 1 1 . = — 1.0. y* — — 2.0. (Case II 5 , or case 5.) 


As an example, one set of ft and y for each (Ii, . . . JIIs) 
is taken and several characteristics corresponding to 
some values of u are given in Figs. 2-16. Table VII 
presents the numerical values of tf x (i— 0, etc. 
Also added is the case E, where ft'= — 0.13, which is, 
of course, included in case III, but this is the actual 
case for the earth. Adopted values of the harmonic 
coefficients are from Kozai (1962b) : 

7 2 = 1082.36X 10- 6 , /*= - 2.566X 10~ 6 , n ~ 
-2.14X10- fi , J h — — 0.063 X 10 -6 . [ } 

Accordingly, 

P= 0.147X10- S fl 0 - 3 , ^=0.0528a 0 - 1 , 

ft — — 0.360 X 10“ 3 <Z(T 2 , n v 

7 = — 2a — 5.540X 10- 3 atr 2 (a= 1 - 5 HVU), 1 ; 

ft"- -4.24 XMHflo- 1 ; 


K 



Fig. 12. 0'*= — 0.25, y* = +2.0. (Case IIIi, or easel.) 


therefore, 

ft' = ft/ ft =-0.129 

7*= — 71.7a 0 2 a-“ 1.99, (7.9) 

h*= 18.9ao/^= 18.9a 0 e sing— 0.0201 

k*= 18.9ao& — 18.9a 0 e cosg, 

where oo stands for the mean semimajor axis, e for the 
eccentricity, and g for the argument of perigee. 

It is noted that in the case of the earth the asym- 
metrical part comes mainly from 7 5 , but not from / 3 . 
The linear shift term (ft") which depends only on / 3 
in our theory is small. It is also noted that ft' is not 
dependent on the mean semimajor axis ao but only on 
/n’s; therefore, if the potential is given, it is an absolute 
constant. 

In concluding this section, we should add some 
remarks on the critical value Uq . As is seen from Table 
VII, the critical value wo corresponds neither to a stable 
point nor to an unstable one except for case 5, where 
of course it corresponds to the unstable points. The 
sign of (gs) o is positive for cases 1-4. Therefore, it might 



Fig. 13. 0'*= -0.25, 7 * = 0.0. (Case III 2 , or case 2" .) 
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be considered that Uo corresponds to a stable point. In 
fact, however, it corresponds to imaginary equilibrium 
points. Or in other words, if we draw the characteristics 
corresponding to m it is split into two circles given by 

7 1 1 

W+[*.TH-W?=--d= , (7.10) 

2 2 (- 2 / 3 ')* 

and for case 5 two circles meet at two points so that 
the two points become unstable points; on the other 
hand, for cases 1-4 they cannot meet each other or the 
radii (us) of one or both circle (s) become (s) imaginary. 



VIII. discussions 

1. As mentioned before in Sec. VI the expansions of 
solutions are so arranged that they are always con- 
vergent. But it is still a problem to find out the most 
efficient expansions for the solutions; in other words, 
to rearrange the series when it is too slowly convergent. 
This situation will occur when q tends to unity. In 
such a case we have another type of expansion by 
exchanging the real period and imaginary period with 
each other, so that for the extreme case of q=l we 
have hyperbolic functions instead of circular functions. 

In any case, in order to have suitable expansions of 
the solutions it is necessary to rearrange them into a 
different form appropriate to each case, respectively. 

Also neglected is a proper method to calculate semi- 
periods (co', a>i', or wi"i) in terms of given constants 



/3', 7 and u. These calculations are closely related with 
the evaluation of q (if A>0) or q x (if A<0). The 
numerical processes to find p and x are also omitted here. 

These practical problems are, of course, important 
if we wish to obtain the solutions in detail but they are 
so complicated that their discussion will be postponed. 

2. The adoption of nondimensional analysis in units 
of the proper powers of 0 is, of course, optional. We 
can treat the calculation without any nondimensional 
analysis; or we may have another type of nondimen- 
sional analysis in terms of 0' or y for example. However, 
the process which we treated in this paper has a slight 
advantage; firstly, in doing so 0 f * is an absolute con- 
stant depending only on the coefficients of harmonics 
but independent of the initial condition; secondly, we 
may treat both the cases /3'>0 or /3'<0, in the same 
framework, for example, we can draw a diagram such 
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Table VII. Critical values of u, and corresponding h* and k*. 


Subcase 

7* 

i 

U*i 

h+i 

k*i 

Remark 

1 

Subcase 

7* 

i 

U*i 

h*% 

k*i 

Remark 




Case I /3*'=— 3.0 






Caselll /?+'= — 0.25 



Ii 

+2.0 

0 

+ 6.2083 

. . . 

. . . 


III! 

+3.0 

0 

+ 0.7656 

. . . 





1 

+ 1.0548 

+0.8847 

0.0 

stable 



1 

+ 0.1371 

+0.2645 

0.0 

stable 

I 2 

+ 1.0 

0 

+ 3.9583 

+1.1072 

0.0 

stable 

III 2 

0.0 

1 


+0.6959 

0.0 

stable 



1 

+ 2 . 0562 

0 

- 0.4844 




2 

+ 0.0734 

-0.8376 

0.0 

stable 








3 

- 0.1295 

-0,2696 

0.0 

unstable 

Ilia 

-1.3 

1 

+ 1.5506 

+1.0111 

0.0 

stable 

I, 

0.0 

1 

0 

2 

+ 3.5139 
+ 2.2083 
+ 1.0702 

+ 1.3008 
-1.1309 

0.0 

0.0 

stable 

stable 



2 

3 

0 

- 0.1717 

- 0.1777 

- 0.2244 

-0.5966 

-0.4145 

0.0 

0.0 

unstable 

stable 



3 

- 0.0840 

-0.1699 

0.0 

unstable 

IH 4 

-1.5 

1 

+ 1.7641 

+ 1.0546 

0.0 

stable 

I* 

-2.0 

1 

+ 7.8553 

+ 1.6290 

0.0 

stable 



2 

- 0.0815 

-0.7296 

0.0 

unstable 



2 

+ 4.6947 

-1.5286 

0.0 

stable 



0 

- 0.1094 



stable 



0 

+ 0.2083 





3 

- 0.1512 

-0.3249 

0.0 



3 

- 0.0500 

-0.1004 

0.0 

unstable 

Ills 

-3.0 

1 

+ 3.9518 

+ 1.3457 

0.0 

stable 

Ib 

-4.0 

1 

+ 14.1383 

+ 1.9056 

0.0 

stable 



2 

+ 1.4070 

-1.1895 

0.0 

stable 



2 

+ 10.3976 

-1.8340 

0.0 

stable 



0 

+ 1.3906 

-1.0000 

±0.6124 

unstable 



0 

+ 0.2083 

-0.0833 

±0.7022 

unstable 



3 

- 0.0775 

-0.1562 

0.0 

stable 



3 

- 0.0358 

-0.0715 

0.0 

stable 




Case E j8+ 

'=-0.13 






Case II /3 

*=-1.0 



Ej 

+3.0 

0 

+ 1.4877 

, . . 

. . . 


Hi 

+2.0 

0 

+ 2.1250 

. . . 




1 

+ 0.0863 

+0.1707 

0.0 

stable 



1 

+ 0.2148 

+0.3855 

0.0 

stable 






II 2 

0.0 


+ 1.0548 

+0.8847 

0.0 

stable 

e 2 

0.0 

1 

+ 0.5269 

+0.6643 

0.0 

stable 

1 

0 

- 0.9573 



0 

+ 0.1250 

. . . 








IIs 

-0.6 

1 

0 

2 

+ 1.6022 

- 0.0850 

- 0.1524 

+ 1.0221 

0.0 

stable 

e 3 

-1.5 

1 

2 

+ 1.6337 
- 0.1397 

+1.0286 

-0.6608 

0.0 

0.0 

stable 

unstable 



-0.6396 

0.0 

stable 



3 

0 

- 0.1655 

-0.3679 

0.0 

stable 



3 

- 0.1698 

-0.3823 

0.0 

unstable . 



— 0.4923 




H 4 

-1.0 

1 

+ 2.0562 

+1.1072 

0.0 

stable 

e 4 

-3.0 

1 

+ 3.7379 

+ 1.3243 

0.0 

stable 



2 

+ 0.0734 

-0.8376 

0.0 

stable 



2 

+ 1.2412 

-1.1618 

0.0 

unstable 



0 

- 0.1250 

. . . 





0 

+ 1.0977 






3 

- 0.1295 

-0.2696 

0.0 

unstable 



3 

- 0.0806 

-0.1625 

0.0 

stable 

II* 

-2.0 

1 

+ 3.5139 

+ 1.3008 

0.0 

stable 

e 6 

-8.0 

1 

+ 18.5553 

+2.0463 

0.0 

stable 



2 

+ 1.0702 

-1.1309 

0.0 

stable 



2 

+ 14.5236 

-1.9847 

0.0 

stable 



0 

+ 0.1250 

-0.2500 

±0.6614 

unstable 



0 

+ 14.5227 

-1.9231 

±0.4865 

unstable 



3 

- 0.0840 

-0.1699 

0.0 

stable 



3 

- 0.0302 

-0.0611 

0.0 

stable 


a 0 hi=ac,e sing 



Fig. 17. /S'* = —0.13, y* = 4-3.0. (Case Ei, or case 1.) 
(Figures 17-21 correspond to the case for the earth.) 

as Fig, 1 on a single sheet; thirdly, the difference 
between 0>O and £<0 can be neglected automatically 
by always selecting the sign of ( p r ] )* or (—/S')* such 


that (£'*)*, (— /3'*)*>0. This is caused by the fact that 
the difference between d>0 and /3<0 remains only on 
the sign of h and k but there is no difference on the 
nature of the critical character, etc. 



Fig. 18. /3'* = — 0.13, 7 *— 0.0. (Case E 2 , or case 2".) 
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SHINKO AOKI 


a 0 hi =a 0 e sin g 



3. We have assumed Z n to be of nth order of magni- 
tude. But if those assumptions are broken down, how 
shall the situation change? We may consider this 
breakdown as follows within our framework : & — > 0, in 
other words a noneffective case of the antisymmetrical 
terms Z 3 and Z 5 (J3 f * — > oo). In this case, the range h* 
and k* for the critical nature (the transition from 
librational to revolutional) is not restricted in a limited 
area according as y* — > — 00 , but h and k themselves 
are limited. Therefore we may let /3 — > 0 without any 
essential difficulty in our theory. 

On the other hand if /S' — > 0, keeping (3 constant the 
range of h* and k * as well as h and k themselves for the 
critical nature are not restricted as 7* — * — oc . There- 
fore, another type of theory would be required which 


a 0 hj =a 0 e sin g 



could treat not only the small eccentricity case but 
also the moderate eccentricity case. In this case where 
the antisymmetrical terms prevail over the sym- 
metrical terms, the general feature, which is constructed 
here by adding the antisymmetrical terms to the 
symmetrical terms or at most by considering the same 
order contributions in some meaning from both sources, 
would be broken down. Fortunately, since 13*' for the 
earth is approximately —0.13, it does not seem to 
necessitate any new theory in this respect. 

4. A higher-order theory beyond the one discussed 
here could be developed along the line which we gave 
in a previous paper (1962). However we must keep in 
mind the restriction made there that the leading terms 
which are, for example, the ones discussed here, should 
prevail over the remaining terms. In other words, if the 
leading terms are so small that the neglected terms 

a o hi=a 0 e sin g 



play an important role in the behavior of the solution, 
then such an approximation process would be broken 
down. In addition, in the present paper, we have taken 
into account two parameters such as (3 and j3\ the 
former is assumed to be third order and the latter is of 
second order; namely the nondimensional quantity 
is assumed to be far from zero or from 

infinity. 

In the previous paper, it is assumed that, for example, 
Z 2 3 , Z2Z4 and Ze are a 11 of the third order of Z2+Z4/Z2; 
this means that, if Z2+Z4/Z2 tends to zero, Z2 3 etc. 
should tend to zero as well, in order to be able to 
apply the general treatment described there, otherwise 
the leading terms are so small that the series obtained 
along the theory would be divergent. A similar situation 
occurs in the present theory in subsection 3. This is 
also true in the case of the higher-order theories here 
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discussed. At any rate it is important to consider first 
the most significant parts of the terms. 

The difference between the previous theory and the 
present one consists of the fact that here the anti- 
symmetrical terms are given the same importance as 
the symmetrical terms by assuming that J n is of the 
wth order of magnitude in the case of small eccentricity. 

IX. CONCLUSION 

In this paper we tried to solve the equations of 
motion of a close satellite near the critical inclination 
in the case of small eccentricity under the influence of 
the potential given by 

J P„(sinS)J. 

Assumed are as follows: the coefficient of zonal har- 
monics J n is of the nth order of magnitude, the eccen- 
tricity is of the first order, and a=l — 5IP/L? is of the 
2nd order. The terms which are taken in the Hamil- 
tonian are of the sixth order of magnitude in this re- 
spect; namely, the terms which have eVa, 2 s / 4 , e?Jz , 
2/2/3, 2/5 etc. as factor are included. It is also noted that 
the lower order terms do not enter the Hamiltonian 
except for unnecessary constant terms. 

In Sec. VI the explicit solutions are given with 
classification according to the interrelation of /3'*, 7* 
and u* [jfor the definition of these quantities, see 
Eqs. (4.4), (2.12) and (3.4)], where jS'* strictly depends 
on the coefficients of zonal harmonics /„, 7* depends on 
the projected angular momentum to the equatorial 
plane as well, and u* is the energy constant in units of 
1 This classification is made in order to have real 


expressions for the solutions. Three parameters involved 
are expressed bv the nondimensional quantities so that 
one may have the actual solutions by factoring a 
proper power of which is assumed as of the first 
order of magnitude. 

For the earth 8* is +0.0528/a 0 ; accordingly, if the 
eccentricity is confined within some range, say 0.05 or 
0.10, then the theory developed here is applicable. 

In any case it is noted that the transition from the 
symmetrical term prevailing case to the antisymmetrical 
prevailing case within the approximation adopted here 
occurs at (between case II and case III, where a great 
change especially in characteristics occurs) : 

This quantity for the earth is approximately —0.13; 
therefore, the earth’s potential is in the anti- 
symmetrical term prevailing case in this respect. 
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